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Cold  collision  shift  cancelation  and  inelastic  scattering  in  a  Yb  optical  lattice  clock 


A.  D.  Ludlow,  N.  D.  Lemke,  J.  A.  Sherman,  and  C.  W.  Oates 
National  Institute  of  Standards  and  Technology,  325  Broadway,  Boulder,  CO  80305,  USA* 

G.  Quemener,  J.  von  Stecher,  and  A.  M.  Rey 

JILA,  National  Institute  of  Standards  and  Technology  and  University  of  Colorado,  Boulder,  CO  80309-0440 

Recently,  p-wave  cold  collisions  were  shown  to  dominate  the  density-dependent  shift  of  the  clock 
transition  frequency  in  a  171  Yb  optical  lattice  clock.  Here  we  demonstrate  that  by  operating  such 
a  system  at  the  proper  excitation  fraction,  the  cold  collision  shift  is  canceled  below  the  5  x  10”ls 
fractional  frequency  level.  We  report  inelastic  two-body  loss  rates  for  3P0  -3Po  and  'S'o  -3Po  scattering. 

We  also  measure  interaction  shifts  in  an  unpolarized  atomic  sample.  Collision  measurements  for  this 
spin-1/2  171  Yb  system  are  relevant  for  high  performance  optical  clocks  as  well  as  strongly-interacting 
systems  for  quantum  information  and  quantum  simulation  applications. 


Large  ensembles  of  ultracold  atoms  offer  atomic  clocks 
a  measurement  of  the  atomic  state  with  high  signal-to- 
noise  ratio.  For  clocks  utilizing  optical  transitions,  this 
has  the  potential  to  yield  time  and  frequency  measure¬ 
ments  with  new  levels  of  precision  and  speed  (e.g.  [1,  2]). 
However,  large  ensembles  of  cold  atoms  can  lead  to  high 
number  density  and  thus  significant  interatomic  interac¬ 
tions.  These  interactions  can  perturb  the  clock  transi¬ 
tion  frequency,  compromising  the  accuracy  of  the  atomic 
standard.  For  example,  in  cesium  fountain  primary  stan¬ 
dards,  cold  collision  shifts  can  become  significant  [3,  4], 
influencing  clock  operation  (e.g.  [5-7]). 

Optical  lattice  clocks,  which  probe  the  ultra-narrow 
'Sq  -3Po  transition  in  two- valence-electron  atoms  held  in 
an  optical  potential,  are  also  susceptible  to  cold  colli¬ 
sions.  Non-negligible  collision  effects  have  been  observed 
in  lattice  clocks  using  nuclear-spin-polarized,  fermionic 
samples  of  87Sr  [8,  9]  and  171  Yb  [10].  Measurement  and 
control  of  these  collisions  therefore  play  a  key  role  in  the 
continued  development  of  these  standards.  At  the  same 
time,  the  control  of  these  interactions  are  an  integral  part 
of  proposals  for  quantum  information  [11,  12]  and  quan¬ 
tum  simulation  of  solid-state-analog  Hamiltonians  [13- 
15].  For  87Sr,  the  primary  interaction  giving  the  cold 
collision  shift  was  identified  as  an  s-wave  interaction  be¬ 
tween  atoms  in  non-identical  superpositions  of  the  clock 
states  [9,  16,  17].  It  was  recently  observed  that  strong 
interactions  can  isolate  and  suppress  the  shift  [16,  18]. 

For  1(1Yb,  it  was  shown  that  a  p- wave  interaction  be¬ 
tween  atoms  in  the  Wo  and  3Po  electronic  states  was  the 
dominant  mechanism  responsible  for  the  cold  collision 
shift  [19].  Unlike  the  s-wave  case,  such  an  interaction  is 
less  sensitive  to  small  particle  distinguishability  between 
the  ultracold  fermions.  Consequently  the  observed  shift 
exhibits  roughly  a  linear  dependence  on  the  3Po  excita¬ 
tion  fraction  (Fig.  1(a)).  Moreover,  in  the  weakly  inter¬ 
acting  regime  of  the  1-D  lattice,  the  cold  collision  shift 
crossed  zero  at  a  mean  excitation  fraction  close  to  0.5. 
In  this  work,  we  exploit  this  zero-crossing  to  demonstrate 
cancelation  of  the  171  Yb  cold  collision  shift  in  a  1-D  opti¬ 


cal  lattice  below  the  5  x  10” 18  fractional  frequency  level. 
The  cancelation  is  enabled  by  the  fact  that  near  50%  ex¬ 
citation,  population  of  the  two  electronic  states  of  one 
atom  induce  equal  collision  shifts  on  the  opposite  states 
of  another  atom,  leaving  zero  net  shift  for  the  clock  tran¬ 
sition  [20].  The  collision  shift  cancelation  can  thus  be 
likened  to  the  Stark  shift  cancelation  which  optical  lattice 
clocks  exploit  by  operating  at  the  “magic”  wavelength 
[21,  22].  Our  uncertainty  in  this  cancelation  reaches  be¬ 
low  the  smallest  total  uncertainty  levels  reported  to  date 
for  any  type  of  atomic  clock,  making  it  a  powerful  and 
pragmatic  technique  for  mitigating  the  collisional  shift  in 
a  lattice  clock. 

The  mYb  (/  =  1/2)  optical  lattice  clock  is  described 
in  detail  elsewhere  [10,  19].  After  dual  laser-cooling 
stages,  atoms  are  loaded  into  a  1-D  optical  lattice  using 
a  laser  at  Amag;c  ~  759  mn.  Atomic  population,  initially 
split  between  the  two  ground  magnetic  substates,  is  opti¬ 
cally  pumped  to  a  single  spin  state  using  er-polarized  light 
resonant  with  the  Wo  (/  =  1/2)  -  3  Pi  (/  =  3/2)  transition 
in  the  presence  of  a  bias  magnetic  field  (5  G).  We  esti¬ 
mate  the  average  atomic  density  can  reach  p0  ss  3  x  1011 
atoms/cm3,  with  an  atomic  temperature  of  ~  10  pK. 
The  atoms  are  then  spectroscopically  probed  on  the  Wo  - 
3P0  transition  using  Ramsey  spectroscopy.  The  Ramsey 
pulse  times  used  here  are  5-10  ms,  while  the  dark  time 
is  80-150  ms.  A  cavity-stabilized  probe  laser  is  actively 
stabilized  to  the  center  Ramsey  fringe  using  an  acousto¬ 
optic  modulator.  The  cold  collision  shift  is  the  measured 
frequency  difference  between  interleaved  atomic  samples 
of  high  and  low  density. 

Measurement  of  the  collision  shift  in  a  1-D  lattice  is 
shown  as  filled  points  in  Fig.  1(a),  as  a  function  of  ex¬ 
citation  fraction  (during  the  Ramsey  dark  time).  The 
theoretical  calculation  of  the  shift  using  a  p-wave  model 
is  also  shown  (solid  line).  In  the  model,  the  domi¬ 
nant  interaction,  Veg,  is  between  Wo  and  3Po  atoms,  and 
a  weaker  interaction  between  two  excited  state  atoms, 
Vee  =  0.1  x  Veg,  is  also  included  [19].  The  shift  has  a 
zero  value  near  50%  excitation,  and  this  zero-crossing  is 
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FIG.  1:  (color  online)  (a)  Cold  collision  shift  as  a  function 
of  excitation  fraction,  for  a  polarized  sample  in  a  1-D  opti¬ 
cal  lattice  oriented  vertically  (filled  triangles)  or  horizontally 
(filled  circles).  Open  squares  are  with  an  unpolarized  atomic 
sample,  (b)  Cold  collision  shift  cancelation:  the  measurement 
uncertainty  for  this  cancelation  is  shown,  given  by  the  total 
deviation  (similar  to  two-sample  Allan  deviation  [23]).  (c) 
Measurement  of  the  residual  cold  collision  shift  when  operat¬ 
ing  at  51%  excitation  and  an  atomic  density  of  pi.  The  red 
dashed  line  indicates  the  weighted  mean  of  the  seven  mea¬ 
surements,  +2.5  mHz,  while  the  solid  lines  gives  the  one- a 
error  bars  of  ±2.4  mHz. 


the  focus  of  our  attention  here.  To  determine  the  precise 
excitation  corresponding  to  zero  shift,  we  made  real-time 
measurements  of  the  excitation  fraction  by  turning  off 
the  second  Ramsey  pulse  and  measuring  the  ground  and 
excited  atomic  populations  after  the  dark  time.  These 
measurements  were  interspersed  between  cycles  of  usual 
Ramsey  spectroscopy.  We  then  slightly  adjusted  the 
probe  laser  power  to  keep  the  measured  excitation  frac¬ 
tion  constant  during  each  shift  measurement. 

The  precision  of  the  zero-crossing  measurement  is  given 
by  the  clock  instability  during  interleaved  measurements 
of  high  and  low  atomic  density.  Our  measurement  sta¬ 
bility  benefits  from  recent  improvements  to  the  cavity- 
stabilized  laser  used  to  probe  the  clock  transition  [1]. 
Fig.  1(b)  shows  the  precision  of  one  such  measurement, 
where  we  observed  a  collision  shift  of  -1.4  (3)  mHz  after 
14500  s  of  averaging.  Fig.  1(c)  shows  the  result  of  seven 
similar  measurements,  taken  sequentially  over  the  course 
of  several  weeks  and  with  different  measurement  dura¬ 
tions.  For  all  but  one  measurement,  the  mean  excitation 
fraction  was  controlled  to  51±0.3±1.3  %,  where  the  first 
uncertainty  is  the  fluctuation  in  the  measured  excitation 
fraction  and  the  second  is  the  systematic  uncertainty  in 


the  absolute  excitation  value.  (For  data  point  number 
two,  the  excitation  fraction  was  set  to  1%  higher;  we 
thus  applied  a  —5  mHz  correction  to  the  measured  col¬ 
lision  shift,  as  determined  from  the  slope  of  the  curve 
in  Fig.  1(a).)  The  weighted  mean  of  the  seven  measure¬ 
ments  (red  lines)  is  2.5  (2.4)  mHz  (reduced  y2  =  1.04). 
This  corresponds  to  a  fractional  shift  of  4.8  (4.6)  xlO-18 
of  the  transition  frequency,  and  demonstrates  the  small¬ 
est  measurement  of  a  collision  shift  in  a  lattice  clock. 

In  order  to  routinely  implement  this  collision  shift  can¬ 
celation,  we  consider  the  robustness  of  this  technique 
to  relevant  experimental  conditions.  For  an  operational 
density  of  pj  ~  3  x  1010  atoms/cm3,  a  1%  change  in  the 
excitation  fraction  leads  to  a  change  in  collision  shift  of 
~  1  x  10-17.  As  described  above,  it  is  straightforward 
to  keep  the  excitation  fixed  at  or  below  the  1%  level. 
However,  a  complication  arises  from  time-dependent  ex¬ 
citation  due  to  trap  loss.  Particularly,  we  have  observed 
inelastic  two-body  losses  involving  both  3Po-3Po  and  (So  - 
3P0  interactions.  With  both  (So  (g)  and  3Po  (e)  popula¬ 
tions  present,  the  number  density  rate  equations  are: 

ng(t)  =  -Tgng(t)  -  (legng{t)ne{t) 
he(t)  =  —&eene(t)2  -  rene(t)  -  pegng(t)ne(t) 

For  a  single  population  ( ng  or  ne)  (3eg  loss  can  be  ignored. 
Fig.  2(a)  shows  ’So  trap  loss  (black  triangles,  for  a  pure 
sample  of  (So  atoms),  with  a  one  body  loss  fit  yielding  a 
trap  lifetime  of  l/rg  =  480  (20)  ms.  Also  shown  is  3Po 
trap  loss  (blue  circles,  for  a  pure  sample  of  3Pq  atoms). 
This  population  experiences  notably  stronger  decay  at 
high  densities,  and  good  fit  requires  both  one-body  (re) 
and  two-body  {(5ec)  losses.  After  integrating  the  losses 
over  the  spatial  extent  of  a  single  lattice  site,  as  well  as 
across  the  distribution  of  occupied  sites,  we  find  1  /Ye  = 
520  (28)  ms  and  pee  =  5  x  10“ 11  cm3/s,  somewhat  larger 
than  the  3Pq-3Pq  decay  measured  in  88Sr  [24,  25].  Because 
the  absolute  atomic  density  is  difficult  to  calibrate,  the 
uncertainty  in  (3ee  is  estimated  as  60%. 

Inelastic  two-body  loss  for  3Po  atoms  can  be  consid¬ 
ered  with  a  single  atom-atom  channel,  time-independent 
quantum  formalism  similar  to  [26,  27].  Lacking  an  accu¬ 
rate  Yb2  potential,  the  short  range  physics  is  described 
here  by  a  boundary  condition  at  an  interatomic  sep¬ 
aration  Ro  =  20  ao  (ao  is  the  Bohr  radius),  repre¬ 
sented  by  two  parameters:  (i)  an  accumulated  phase 
shift,  0  <  5  <  it,  due  to  the  unknown  short-range  po¬ 
tential  from  R  =  0  to  R  =  Rq,  and  (ii)  a  loss  probability, 
0  <  pis  <  1  at  R  =  I?o-  The  long  range  physics  is  ac¬ 
curately  described  by  a  van  der  Waals  interaction,  with 
Ce  =  3886  a.u.  given  by  [28],  The  Schrodinger  equa¬ 
tion  is  numerically  solved  from  Rq  to  R  — >  oo,  and  the 
cross  section  is  computed  as  a  function  of  the  collision 
energy.  The  thermalized  rate  coefficients  are  found  by 
averaging  the  cross  sections  with  a  Maxwell-Boltzmann 
distribution  of  collision  energies  for  a  given  temperature. 
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FIG.  2:  (color  online)  (a)  'So  trap  loss  (black  triangles)  with 
an  exponential  fit  (dashed  curve).  3P0  trap  loss  (blue  circles) 
with  a  one-  and  two-body  loss  fit  (solid  curve).  The  inset 
shows  'So  loss  for  a  pure  'So  sample  (black  triangles)  and  a 
mixed  sample  of  'So  and  3Po  atoms  (blue  circles).  The  former 
is  fit  with  a  simple  exponential,  the  latter  with  3Po  -3Po  and 
'So  -3P0  two  body  losses,  (b)  Two  body  loss  rates  as  a  function 
of  temperature,  calculated  with  p\s  =  1  (dashed  curves)  and 
with  pis  =  0.8  and  S  =  0.51  n  (solid  curves).  Black  is  for 
a  nuclear-spin-polarized  sample  of  3Po  atoms  (See),  red  an 
unpolarized  sample  of  3Po  atoms  ([See  +  0ee\/2).  Circles  give 
experimental  measurements,  black  for  See  and  blue  for  /3eg. 

With  full  loss  probability  (p^  =  1),  /3ee  is  shown  in  Fig. 
2(b)  (black  dashed  curve)  as  a  function  of  temperature. 
Here,  the  rates  are  universal,  independent  of  <5  [26,  27]. 
The  interaction  involves  two  identical  fermions,  and  thus 
quantum  statistics  dictates  that  the  interactions  must  be 
odd  partial  waves,  notably  p-wave  at  these  ultracold  tem¬ 
peratures.  As  two  excited  atoms  approach  each  other, 
those  which  successfully  tunnel  through  the  p-wave  bar¬ 
rier  are  assumed  to  inelastically  scatter  with  unit  prob¬ 
ability,  a  reasonable  assumption  due  to  the  large  num¬ 
ber  of  exit  collision  channels  available  to  a  pair  of  atoms 
with  high  internal  energy.  For  an  unpolarized  sample 
of  3P0  atoms  (here  with  equal  my  =  itl/2  populations), 
interactions  between  distinguishable  atoms  also  include 
an  s-wave  term,  and  the  loss  rate  for  distinguishable  3Po 
atoms  is  labeled  $ee.  The  total  loss  rate  for  the  unpo¬ 
larized  sample,  given  by  the  average  of  f5ee  and  /3ee,  is 
shown  in  Fig.  2(b)  as  a  red  dashed  curve  (p\s  =  1).  Ex¬ 
perimentally,  we  observed  identical  loss  rates  within  10% 
for  a  polarized  and  an  unpolarized  sample  at  10  pK  (Fig. 
2(b),  black  circle).  Since  full  loss  predicts  unequal  loss 
rates  for  polarized  and  unpolarized  samples  at  10  pK,  3Pq 
atoms  may  not  be  lost  with  full  unit  probability  at  short 
range.  Instead,  using  short  range  parameters  of  pis  =  0.8 


and  6  =  0.51 7r,  the  loss  rates  are  shown  as  solid  curves 
in  Fig.  2(b)  for  polarized  (black)  and  unpolarized  (red) 
cases.  Here,  we  found  better  agreement  with  the  exper¬ 
imental  data,  suggesting  a  deviation  from  the  universal 
regime. 

For  a  mixed  population  of  1Sq  and  3Po  atoms,  we  ob¬ 
served  additional  loss  through  f3eg.  Fig.  2(a)  inset  high¬ 
lights  this  by  comparing  '.So  trap  loss  for  two  different 
atomic  samples:  with  (blue  circles)  and  without  (black 
triangles)  the  presence  of  3Po  atoms.  The  additional  loss 
is  particularly  notable  at  short  times  where  both  '.So  and 
3P0  populations  are  large  and  is  consistent  with  a  'S’o  -3Po 
two-body  loss  rate  at  the  level  of  f3eg  =  3  x  10_n  cm3/s. 
The  3P0  population  is  excited  coherently  from  the  '.S’o 
state,  but  due  to  excitation  inhomogeneity,  the  lossy  col¬ 
lisions  are  between  partially  distinguishable  atoms  [19]. 
The  magnitude  of  the  inelastic  loss  is  noteworthy  because 
for  molecular  states  correlating  to  a  ^Sq  and  3Po  atom 
pair,  only  the  *E+  ground  state  (correlating  to  the  'So- 
'So  state)  lies  at  lower  energy,  and  long-range  coupling 
to  this  state  is  spin-forbidden. 

In  general,  all  of  these  loss  processes  lead  to  a  time- 
dependent  excitation  fraction  during  the  Ramsey  dark 
time.  This,  in  turn,  affects  the  balance  between  'So  and 
3P0  collisionally  induced  energy  shifts.  For  simplicity,  we 
can  easily  operate  at  a  lower  atom  number  density,  reduc¬ 
ing  not  only  the  collisionally-induced  shifts,  but  also  the 
two-body  inelastic  losses.  At  the  same  time,  the  number 
of  quantum  absorbers  remains  well  in  excess  of  1000,  in 
order  to  accommodate  a  high  signal-to-noise  ratio.  At 
a  density  of  pi,  the  excitation  fraction  over  a  dark  time 
of  T  =  150  ms  changes  by  only  several  percent.  We 
therefore  take  the  time-averaged  value  to  indicate  the 
excitation  fraction  where  the  shift  is  canceled. 

Any  degree  of  imperfect  polarization  of  the  nuclear 
spin  state  introduces  a  host  of  other  possible  atomic  inter¬ 
actions.  In  p-wave,  the  V~  interaction,  which  shifts  the 
singlet  state,  becomes  allowed  [19].  Furthermore,  distin- 
guishability  between  different  nuclear  spin  states  allows 
s-wave  interactions  Ugg,  Uee,  and  17+ .  In  1  ‘ 1  Yb,  both  the 
s-  and  p-wave  |gp)  interaction  terms  are  small  [29].  How¬ 
ever,  the  remaining  interactions  can  alter  the  observed 
collision  shift  from  the  spin-polarized  case.  For  this  rea¬ 
son,  high  polarization  purity  is  important  for  optical  lat¬ 
tice  clocks.  Here,  we  benefit  from  the  simple  structure 
( I  =  1/2)  of  mYb  and  can  readily  optically  pump  to  a 
single  spin  ground  state  with  99%  purity.  The  polariza¬ 
tion  purity  is  directly  measured  by  observing  the  absence 
of  a  clock  excitation  spectrum  for  the  unpopulated  spin 
state. 

To  quantify  the  effect  of  imperfect  polarization,  we 
measure  the  collision  shift  as  a  function  of  excitation  for 
unpolarized  atoms  (open  squares  in  Fig.  1(a),  equal  mix¬ 
ture  of  both  spin  states).  Since,  during  spectroscopy, 
a  weak  bias  field  (1  G)  lifts  the  degeneracy  of  the  n- 
transitions  from  the  two  spin  states,  here  the  excitation 
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FIG.  3:  (color  online)  Cold  collision  shift  as  a  function  of 
excitation  for  a  polarized  sample  in  a  2-D  lattice,  for  different 
Ramsey  times.  p2  is  a  typical  2-D  lattice  density,  where  25% 
of  the  atoms  are  doubly  occupied  in  a  lattice  site  with  density 
of  4  x  1012  atoms/cm3.  Circles  correspond  to  experimental 
measurements,  and  solid  curves  are  theoretical  calculations 
[19].  In  both  cases,  blue  is  a  Ramsey  time  of  T  =  10  ms.  red 
is  40  ms,  black  is  80  ms,  green  is  160  ms,  and  orange  is  210 
ms. 


fraction  is  defined  as  that  for  the  particular  spin  state 
being  resonantly  excited,  not  accounting  for  the  other 
‘spectator’  spin  state.  In  general,  the  measured  shifts  are 
smaller  than  for  the  spin-polarized  case,  implying  that 
competing  interactions  have  the  opposite  sign  as  those  in 
the  polarized  case.  Notably,  the  measured  zero-crossing 
in  the  shift  occurs  at  a  lower  excitation  fraction  (around 
40%),  leaving  a  net  positive  shift  at  51%,  where  the  shift 
is  zero  for  the  spin  polarized  case.  Based  upon  these 
measurements,  we  determine  a  1%  polarization  impurity 
will  not  affect  the  shift  zero-crossing  at  or  above  the  level 
measured  in  Fig.  1(c). 

We  have  shown  that  the  cold  collision  shift  can  be  can¬ 
celed  at  the  5  X 10“ 18  level  in  a  1-D  optical  lattice  clock  of 
mYb.  Other  implementations  may  be  suitable  for  reduc¬ 
ing  the  clock  shift  at  or  below  this  level.  In  a  2-D  optical 
lattice  clock,  strong  interactions  can  lead  to  a  decay  of 
the  collision  shift  [19,  30].  In  particular,  longer  Ramsey 
dark  times  lead  to  smaller  shifts  and  can  reduce  the  shift 
dependence  on  excitation  fraction  (Fig.  3),  making  it  at¬ 
tractive  for  shift  reduction.  Care  must  be  taken  since,  as 
we  have  observed  both  experimentally  and  theoretically, 
the  interactions  can  also  reduce  Ramsey  fringe  contrast. 
Alternatively,  the  2-D  lattice  system  also  exhibits  a  zero 
crossing  in  the  shift  versus  excitation.  However,  strong 
interactions  can  move  the  zero-crossing  excitation  with 
a  nonlinear  dependence  on  interaction  strength  (Fig.  3), 
and  thus  the  weak  interactions  of  the  1-D  lattice  may  be 
easier  to  control. 

The  3-D  optical  lattice  continues  to  be  an  interesting 
choice  [21],  where  it  is  straightforward  to  achieve  high 
atom  number  with  an  average  of  <  1  atom  per  lattice 
site.  The  small  fraction  of  lattice  sites  with  double  oc¬ 
cupancy  will  exhibit  very  strong  interactions.  Doubly- 
occupied  sites  could  also  be  eliminated  using  photoas¬ 


sociation  losses  [31]  or  directly  via  the  two-body  losses 
observed  here.  Furthermore,  the  3-D  optical  lattice  sys¬ 
tem  may  exhibit  kinematic  suppression  of  the  atomic  in¬ 
teractions  which  yield  the  shifts  [32].  Indeed,  these  rich 
atomic  systems  will  undoubtedly  continue  to  offer  many 
interesting  phenomena  in  1-,  2-,  and  3-D  confinement. 
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